Journal of Applied Crystallography

Accurate determination of state occupancies is essential for interpreting the structural heterogeneity inherent in time-resolved crystallography. However, in cases of high spatial overlap between states, as commonly observed in time-resolved crystallography data, the strong correlation between occupancy and atomic displacement parameters (ADPs) can render single point estimates from standard refinement protocols unreliable. We introduce MEROS (Multi-state Ensemble Refinement for Occupancy Statistics), a pipeline that implements an ensemble refinement approach to assess the post-refinement occupancy-ADP statistics of multiple overlapping states. MEROS utilizes a Monte Carlo sampling of the parameter space, performing independent refinements from randomized starting occupancies and ADP values to empirically characterize the convergence and uncertainty of the solution. The method is implemented as a modular Python pipeline that wraps established refinement programs, ensuring compatibility with existing workflows. We demonstrate its applicability in two case studies: a two-state ligand binding model in T4 lysozyme L99A and a four-state covalent catalysis mechanism in {beta}-lactamase CTX-M-14. MEROS provides occupancy and ADP mean values with standard deviations that directly quantify the informational content of the experimental diffraction data.

Due to recent technological advances, in situ structural cell biology is becoming a high throughput microscopy technique as all the steps of the workflow, from sample preparation to data analysis, are executed faster, more reliable and more reproducible. Sample thinning by cryoFIB-SEM is an essential tool in preparing electron transparent lamellae of biological specimens suitable for further characterization by cryoET. Modern cryoFIB-SEM instruments can be operated remotely and are capable of automated and unsupervised lamellae preparation. To take full advantage of these developments they need a constant supply of LN2 to maintain cryogenic conditions inside the microscope chamber. Here, we introduce a custom automated LN2 refill system that is compatible with gas cooled cryostages, supports long-term cryoFIB-SEM operations and liberates the user from highly repetitive and manual work. We believe this solution can be utilized with other cryoSEM or cryoFIB-SEM devices requiring N2 gas-flow cooling.

Structure elucidation of biological macromolecules by single particle cryogenic electron microscopy (SPA cryo-EM) or cryogenic electron tomography (cryo-ET) relies on low-dose imaging on cryogenic transmission electron microscopes (cryo-TEMs). Routine microscope setup remains technically demanding and can be time-consuming, particularly for inexperienced or infrequent users. We present LowDoseWizard, a guided workflow implemented in SerialEM that enables rapid and standardised setup of cryo-TEM imaging conditions. From minimal user input, the workflow configures microscope optics, camera parameters and image shift settings for all low-dose imaging states, and guides the user through key daily alignment procedures including beam shift offset calibration, objective lens astigmatism correction and coma-free alignment. The workflow is organised into modular routines that can be executed sequentially or independently, while microscope-specific acquisition parameters are defined in editable configuration files, allowing flexible adaptation to different instruments without modification of the core scripts. Across user sessions on three microscopes at EMBL Heidelberg, the complete setup required on average less than 15 minutes. To assess whether predefined imaging conditions generated by the workflow are compatible with high-resolution data collection, we acquired apoferritin data on a 200 kV Glacios and a 300 kV Titan Krios. These datasets yielded reconstructions at 1.62 [A] and 1.09 [A] resolution, respectively, demonstrating that rapid, guided setup can support near-atomic and atomic-resolution single particle cryo-EM. LowDoseWizard lowers the barrier to robust cryo-TEM setup, reduces the time spent on routine parameter selection and alignment, and helps users focus on sample-specific aspects of data acquisition such as target selection. The workflow should be particularly valuable in shared instrumentation environments, where accessibility, reproducibility and efficient microscope use are critical.

Quantum kernel methods offer a potential advantage for classification tasks in high-dimensional feature spaces, yet their practical benefit critically depends on how input features are prepared. We compare five dimensionality reduction strategies--principal component analysis (PCA), Gaussian random projection (RP Gaussian), sparse random projection (RP Sparse), partial least squares (PLS), and uniform manifold approximation and projection (UMAP) -- as pre-processing steps for quantum kernel support vector machines (SVMs) applied to trabecular bone classification from synthetic micro-computed tomography (micro-CT) data. Using a custom procedural generator based on Gaussian random field zero-crossings, we produced 500 synthetic trabecular bone volumes with controlled morphometric properties such as bone volume fraction (BV/TV), trabecular thickness (Tb.Th), number (Tb.N) and spacing (Tb.Sp). Texture features extracted from grayscale slices are reduced to 8-dimensional quantum circuit inputs via each method, then classified using both classical radial basis function (RBF)-SVMs and quantum kernel SVMs with ZZ feature maps on a statevector simulator, both evaluated with 5 x 5 repeated stratified cross-validation (25 folds). Our results show that UMAP is the only reduction method where the quantum kernel remains competitive with the classical baseline. Under repeated cross-validation, UMAP showed a +0.032 accuracy gap favouring the quantum kernel (Dietterich 5 x 2 CV p = 0.177); however, validation on 10 fully independent datasets--each with independently generated samples, separate reduction fits, and separate kernel matrices -- reversed the sign to -0.030 (paired t-test p = 0.123; Wilcoxon p = 0.193; quantum wins 3/10 datasets), indicating that the apparent advantage was likely inflated by fold dependence. Nevertheless, UMAPs gap remains small and non-significant in both analyses, whereas all linear methods (PCA, RP Gaussian, PLS) show substantial quantum deficits of -0.090 to -0.116 across BV/TV classification, with PCA and PLS remaining significant under corrected tests (5 x 2 CV p = 0.004 and p = 0.007 respectively). We additionally evaluate quantum kernel ridge regression for continuous morphometric prediction, finding that ZZ quantum kernels fail uniformly at regression (negative R2 for all methods except PLS at 4 qubits), suggesting that the ZZ kernel captures decision boundaries but not smooth metric structure. These findings provide practical guidance for feature engineering in near-term quantum machine learning pipelines and demonstrate that the choice of dimensionality reduction can determine whether quantum kernels remain competitive with classical baselines.

GFAP is a type III intermediate filament primarily found within astrocytes and is known to maintain proper cell structure and mechanical strength. Mutations in GFAP are implicated in the pathology of Alexander disease, a neurodegenerative disease characterized by cytoplasmic inclusions of protein, known as Rosenthal fibers. GFAP has a typical type III intermediate filament domain structure, consisting of a highly conserved alpha-helical rod domain bracketed by an intrinsically disordered N-terminal head and C-terminal tail domains. While the general domain organization of monomeric GFAP and the assembly process for higher order quaternary structures are known, we lack an atomic resolution mechanistic understanding of GFAP assembly into mature filaments. Understanding the structure of GFAP filaments and how mutations disrupt this structure will provide vital information into how mutations produce Alexander disease pathology. As a first step towards a mechanistic description, we characterized GFAP wild type tetrameric and filamentous assemblies using solid state NMR and compared the results to those obtained from an assembly-deficient GFAP mutant. For wild-type GFAP, we observe surprisingly uniform rigid alpha helical structure and can spectroscopically resolve highly mobile intrinsically disordered regions in the filament assemblies. Wild type tetramers show increased mobility, likely arising from the head and tail domains. Mutation of the highly conserved cysteine at position 294 to serine results in an inability to form full-length filament assemblies. We show that the rigid regions of the C294S mutant assemblies largely remain structurally consistent with wild type tetrameric assemblies but differ from wild-type filament assemblies. There is an increase in highly mobile regions for the C294S mutant relative to the wild-type. Our results provide a foundation for developing solid state NMR approaches to characterize intermediate filament assembly mechanisms and the interfering effect of disease mutations.

Aims/hypothesisQuantitative nanoscale analysis of insulin secretory granules (ISGs) in human pancreatic tissue has been limited by the lack of imaging methods that combine high resolution with large-scale sampling. We aimed to establish expansion microscopy (ExM) as a platform for in situ, quantitative analysis of ISG organisation in human {beta}-cells and to assess whether type 2 diabetes (T2D) is associated with alterations in granule size, abundance or spatial organisation. MethodsWe applied Magnify ExM to PFA-fixed, paraffin-embedded pancreatic tissue sections from 6 human donors, 3 non-diabetic (ND) and 3 T2D, enabling super-resolution optical imaging of insulin-labelled granules. Insulin-positive structures were segmented and analysed using a morphometric pipeline to quantitatively assess size, shape and spatial features. Granule clustering was quantified based on combined area and roundness criteria. ResultsThe diameter distribution of highly circular granules was similar between ND and T2D samples and estimates of granule number per cell indicated only a modest reduction in T2D ([~]25%). In contrast, mapping insulin-positive structures in a roundness-area space revealed a marked enrichment of large, irregular objects consistent with granule clustering in T2D. The fraction of clustered granules was significantly increased in T2D and strongly inversely correlated with insulin stimulation index (r = -0.85). Conclusions/interpretationThese results establish expansion microscopy as a powerful platform for quantitative nanoscale analysis of human pancreatic tissue and identify altered spatial organisation of insulin granules, rather than marked granule depletion, as a prominent feature associated with {beta}-cell dysfunction in T2D. Research in contextO_ST_ABSWhat is already known about this subject?C_ST_ABSO_LI{beta}-cell dysfunction in type 2 diabetes is often attributed to reduced insulin content or {beta}-cell loss. C_LIO_LIInsulin secretory granules (ISGs) have been characterised ultrastructurally, but quantitative analysis in human tissue remains limited. C_LIO_LISuper-resolution approaches, including expansion microscopy, are emerging tools for nanoscale imaging in biological tissues. C_LI What is the key question?O_LIIs {beta}-cell dysfunction in type 2 diabetes associated with depletion of insulin granules or with altered spatial organisation? C_LI What are the new findings?O_LIInsulin granule size distribution is largely preserved in type 2 diabetes, with only a modest reduction in granule number per cell. C_LIO_LIA significant increase in insulin granule clustering is observed in diabetic {beta}-cells. C_LIO_LIGranule clustering is strongly inversely correlated with insulin secretion in the same donor tissues. C_LI How might this impact on clinical practice in the foreseeable future?O_LIIdentifying altered granule organisation as a feature of {beta}-cell dysfunction may help refine the understanding of disease mechanisms and guide future strategies targeting {beta}-cell function. C_LI

Ionizing radiation induces molecular responses that may be used to estimate exposure when physical dosimeters are unavailable. Here we present two large-scale proteomics datasets generated from mouse dorsal skin punch samples collected following controlled X-ray exposures spanning multiple doses, dose rates, and post-exposure time points. Experiment 1 comprised 96 samples (including 16 reference samples) collected 6 days after exposure to 0-75 cGy delivered at either 30 or 300 cGy/min. Experiment 2 comprised 936 samples (including 236 reference samples) exposed to 0-100 cGy at either 3 or 28 cGy/min dose rates and harvested between 7 and 150 days post-exposure. All samples were processed using a standardized workflow involving automated bead-based digestion and data-independent acquisition mass spectrometry. The datasets include multiple pooled reference sample types, process controls, and system suitability standards ensuring high quality data. All data presented are available via ProteomeXchange at several levels of processing, from raw files through normalized peptide- and protein-level abundance matrices suitable for biomarker discovery and machine learning applications. This dataset will facilitate generation of new insights into the biological changes and molecular signatures resulting from X-ray exposure in mice and may also help inform future studies in humans.

Invertebrate vision relies on bistable visual pigments flipping upon photon absorption between rhodopsin and metarhodopsin states. In living butterflies, the UV-VIS absorption spectra of rhodopsin and metarhodopsin, respectively with 11-cis and all-trans isomers of 3-hydroxy-retinal (A3) chromophore, can be conveniently recorded from the eyeshine, the light reflected from the compound eye after passing twice through the light-guiding rhabdoms. * Here, a microscope coupled with a broadband LED source and a microspectrometer was used to record photorelaxations reported in eyeshine reflection spectra. Fitting temporal exponential relaxations to log-reflectance arrays yielded transient and baseline spectra that are analogous to absorbance difference and sum, respectively. Both types of spectra were subjected to singular value decomposition and to fitting of templated visual pigment absorption spectra. * The compound eye of the high brown fritillary Fabriciana adippe was exposed to a series of second-long broadband light pulses, causing photorelaxations with time constants between 40 and 120 ms that led to 80% metarhodopsin in equilibrium. The transient and baseline spectra were fitted with pigment templates, estimating the alpha peak wavelength 547-552 nm for rhodopsin and 496-501 nm for metarhodopsin. The metarhodopsin to rhodopsin alpha peak absorbance ratio 1.25-1.35 is consistent with the isosbestic wavelength at 530 nm. The second isosbestic wavelength indicates that rhodopsin beta (UV) peak absorbs more strongly than metarhodopsin below 405 nm. * Baseline spectra, which were not explicitly analysed in previous studies, enable concatenation of exposures, monitor long-term changes of pigment, and enhance the estimation of beta peak parameters. * The method can be directly used in many butterflies and could be adapted to other insects, particularly fruitflies, facilitating studies of the relation between the visual pigment spectra and the opsin sequences. Spectroscopic results can be complemented with physiologically measured photoreceptor spectral sensitivity datasets and analysed with the same global fitting procedure.

Species of Bacillus bacteria including Bacillus safensis and Bacillus subtilis are finding increasing uses in biotechnology and bioremediation, thanks in part to their metabolic robustness. Malate dehydrogenase (MDH) is at the heart of central metabolism and thus a better understanding of Bacillus MDH proteins could aid in the optimization of these applications. MDH of Bacillus spp. belong to the lactate dehydrogenase (LDH)-like class of MDHs, otherwise known as the MDH3 class. Despite wide prevalence in nature among prokaryotes and archaea, this typically homotetrameric class is understudied compared to the MDH1 and MDH2 classes found in eukaryotes. We therefore recombinantly expressed and purified MDH proteins from two societally relevant Bacillus spp.-B. safensis and B. subtilis-and characterized them biophysically (via Size Exclusion Chromatography-Small Angle X-ray Scattering (SEC-SAXS) and Differential Scanning Fluorimetry (DSF)) and enzymatically (via spectroscopic activity assays). As expected based on their high sequence identity, the two MDH orthologs had similar properties in most regards, including a tetrameric structure and high susceptibility to substrate inhibition. However, we uncovered differences in conditional thermal stability, in addition to subtle differences in enzymatic activity that offer insight into the workings of LDH-like MDH. Summary statementMalate dehydrogenase (MDH) is a fundamental metabolic enzyme, from microbes to mammals, yet comparably little is known about microbial MDH, especially MDH of the tetrameric MDH3 class. We compare the biophysical and enzymatic properties of two such enzymes from the societally relevant bacterial species Bacillus subtilis and Bacillus safensis, offering useful insight with potential biotechnological implications.

We introduce a workflow to identify oligomeric structures that are recorded with single-molecule localization microscopy (SMLM) under cryogenic conditions. Typically, these oligomers are assumed to consist of protomers arranged as equilateral two-dimensional polygons and every protomer is labeled with a dye molecule for visualization. Unlike previous work, we consider scenarios in which the sample plane has an unknown orientation relative to the focal plane. Our contribution is a high-precision plane-fitting algorithm to determine the sample plane, combined with geometrical transformations and two circle-fitting algorithms to identify the oligomeric structures. Our simulations on synthetic data demonstrate that the proposed workflow achieves high accuracy in estimating both the unknown tilted plane and the oligomer size.

Stimulated emission depletion (STED) microscopy is a super-resolution fluorescence imaging technique that achieves high spatial and temporal resolution by exploiting stimulated emission to induce fluorescence depletion (FD) and is expected to have substantial utility for imaging applications using fluorescent proteins. However, the compatibility of fluorescent proteins with STED microscopy systems has been understood primarily through empirical observations, and there is no established methodology for the rational selection of fluorescent proteins for STED microscopy. In this study, we systematically evaluated the compatibility of commonly used fluorescent proteins with STED microscopy systems by measuring FD properties using transient absorption spectroscopy and fluorescence dip spectroscopy, both of which are classified as two-color spectroscopy (TCS). Fluorescent proteins identified as compatible with the STED microscopy system based on the TCS measurements were employed for three-dimensional STED imaging of cellular samples expressing each protein. In all samples, three-dimensional spatial resolution was improved relative to confocal laser microscopy, with particularly marked improvements in z-axis resolution. These findings demonstrate that measurements of FD properties via TCS provide a robust approach for evaluating the compatibility of fluorescent proteins with the STED microscopy system and for selecting suitable fluorescent proteins for STED imaging.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, can use host derived lipids as carbon and energy source for survival. Mammalian cell entry (Mce) associated membrane (Mam) proteins are important for the stability of lipid importing Mce complexes. Mtb has five homologs of Mam proteins referred as orphaned Mam (OmamA-E) proteins. A recent study suggested that OmamC (Rv1363c) is essential for the storage and utilization of lipids under starvation in Mtb. To understand the structure and interactions of OmamC, we generated a truncated soluble variant of OmamC (OmamC129-261). Here, we report on the challenges encountered during the crystallization and structure determination of OmamC129-261 and the strategies applied to overcome them. Despite the AlphaFold2 predicted model proving an initial molecular replacement solution, experimental phasing was necessary to determine the structure of OmamC129-261. Heat treatment of protein prior to crystallization setup removed partially unfolded protein present and played a critical role in enhancing the reproducibility and diffraction quality of OmamC129-261 crystals. Although reported earlier, it is not a widely used method. It is worth to try this method, especially, when faced with poor reproducibility and diffraction of crystals.

Human milk contains structurally diverse glycans with key roles in shaping infant development, yet analytical constraints limit characterisation from low-volume samples. Glycosaminoglycans (GAGs), including chondroitin sulphate (CS), are understudied due to existing protocols requiring sample volumes of at least 5 mL and lengthy extraction steps prior to instrumental analysis. This study establishes a workflow for quantifying CS disaccharides from 25 {micro}L of human milk, enabling analysis of samples previously inaccessible to GAG profiling, such as those collected as salvage samples from neonatal intensive care units. For CS quantification, the CS is first enzymatically depolymerised using chondroitinase ABC to release repeating disaccharide units. Matrix complexity is reduced via two rounds of acetonitrile-based protein and lipid precipitation. Disaccharides are separated by hydrophilic interaction liquid chromatography and detected using a Triple Quadrupole Mass Spectrometer, providing robust sensitivity for all CS disaccharides. Method development and validation were performed using pooled mature human milk from term infants. This workflow facilitates detection of all CS disaccharides, with low but reproducible recoveries for total CS. Low- and high-level spike recoveries were 41.3% (RSDr 7.5%, RSDiR 15.9%) and 43.7% (RSDr 24.4%, RSDiR 27.9%), respectively. Despite modest absolute accuracy, precision remained sufficient to make relative comparison of CS concentrations between samples. This method expands the analytical toolkit for human milk glycomics, enabling same day preparation and CS profiling from sample volumes that are 200 times smaller than prior work, supporting future investigations into GAG-mediated functions in early life. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=134 SRC="FIGDIR/small/723732v1_ufig1.gif" ALT="Figure 1"> View larger version (31K): org.highwire.dtl.DTLVardef@176dffborg.highwire.dtl.DTLVardef@16ae4ccorg.highwire.dtl.DTLVardef@d333c2org.highwire.dtl.DTLVardef@1eb3216_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical abstractC_FLOATNO Schematic of sample preparation protocol 25 L of human milk is combined with lyase enzymes and TRIS buffer containing the internal standard prior to incubation. Samples then undergo multiple rounds of centrifugation and refrigeration before analysis via LC-MS/MS. Made using BioRender.com. Glycan nomenclature following Varki et al., 2015. C_FIG

Coarse-grained structure-based models (CG-SBMs; or G[o] models) are simplified potential energy functions of biomolecules or biomolecular complexes that encode their structure. Molecular dynamics simulations of such SBMs have been successfully used to study long time-scale dynamics such as protein and RNA folding, and large conformational transitions of biomolecular complexes. SBMs have several advantages: (1) Their MD simulations are computationally inexpensive, making extensive sampling easily accessible to many researchers. (2) They are easy to modify and can be adapted for the specific biomolecular problem that needs to be investigated. However, the force-fields of SBMs are not usually included in commonly used biomolecular simulation packages resulting in a barrier to their use. Here, we present SuBMIT (Structure Based Models Input Toolkit; https://github.com/sglabncbs/submit), a toolkit for generating coarse-grained SBM input files for performing MD simulations with GROMACS and OpenMM/OpenSMOG. Simulations whose input files can be generated using the different flavors of CG-SBMs present in SuBMIT include the folding and conformational ensembles of proteins with intrinsically disordered regions, 3D-domain-swapping in proteins and the dynamics of RNA-protein assemblies (e.g., simple RNA viruses).

Developing parsimonious, mechanism-aware quantitative models that predict how biological effectiveness changes with different modifiers remains, in general, an unsolved problem. Advances in radiobiological research have created a large knowledge base of first-principles mechanistic models of radiation response that, in principle, could accurately predict radiosensitivity across different experimental and clinical conditions. However, in practice these mechanistic models come with an overabundance of parameters, the majority of which are practically unidentifiable and, moreover, likely unnecessary if one simply wishes to predict how radiosensitivity changes for some specific modifier of interest. Nevertheless, determining which few details in the full mechanistic model are relevant for a given purpose, as well as how to remove any other extraneous details, remains a highly non-trivial task. In this study, we demonstrate the potential of model reduction, starting from a detailed mechanistic description, as a systematic strategy for deriving parsimonious, experimentally falsifiable radiobiological descriptors. As a proof-of-concept demonstration, we apply the Manifold Boundary Approximation Method (MBAM) to a Mechanistic Model of DNA Repair and Survival (MEDRAS), for the problem of cell survival prediction following an acute exposure. Our findings reveal that the complete MEDRAS model for an arbitrary mixed-quality exposure can be structurally simplified to a reduced three-parameter model for an effective uniform-quality, named MEDRAS-LPL. Additional MBAM analysis on MEDRAS-LPL identifies two boundaries in parameter space, corresponding to sparsely ionizing and densely ionizing radiation. Mapping of MEDRAS-LPL parameter space on to effective LQ space further demonstrates that parameters close to the sparsely ionizing boundary line up with expectations from the theory of dual radiation, while parameters close to the densely ionizing boundary line up with expectations from a purely linear model based on a target-theory description. Moreover, our formalism predicts enhanced synergistic interactions between sparsely ionizing and densely ionizing radiation beyond the Zaider Rossi model (ZRM) paradigm, in line with empirical observations. The results highlight the potential for using reduced-order models not only for predictive applications but also for generating novel hypotheses that can inform future experimental designs and optimization strategies in radiobiology.

The Tar-DNA Binding Protein-43 C-terminal region, TDP43LC, has been previously shown to form amyloid-like fibrils with distinct folds in ALS and FTD. In both diseases, proteinaceous inclusions contain TDP43 C-terminal protein fragments as well as phosphorylated TDP43. Here, we use solution NMR to show that soluble phosphomimetic TDP43LC, P-TDP43LC, is structurally similar to wild-type TDP43LC. Disperse P-TDP43LC, like wild-type protein, contains a central helical region flanked by long disordered regions. Despite this similarity, our turbidity measurements, imaging, and kinetic assays show that P-TDP43LC has different aggregation behavior than wild-type protein. Using solid state NMR measurements we find that that phosphomimetic mutations alter the wild-type fibril conformation. Electrostatic repulsion from negatively charged sidechains, despite having little effect on the soluble proteins structure, perturbs amyloid-like fibril formation and selects for a different conformation in vitro. These results shed light on the structural role of TDP43LC phosphorylation in fibril formation in disease. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=104 SRC="FIGDIR/small/725298v1_ufig1.gif" ALT="Figure 1"> View larger version (16K): org.highwire.dtl.DTLVardef@1c63aforg.highwire.dtl.DTLVardef@1d48ed6org.highwire.dtl.DTLVardef@1ed8fd3org.highwire.dtl.DTLVardef@17d67a8_HPS_FORMAT_FIGEXP M_FIG C_FIG SynopsisPhosphomimetic mutations at ALS and FTD neurodegeneration-associated sites in an amyloid forming protein perturbs the aggregated structure compared to wild-type protein.

Metabolomics offers a sophisticated analytical framework for characterising the molecular phenotype of biological organisms and complex living systems at a high resolution. As the functional endpoint of the omics cascade, the metabolome serves as a close reflection of cellular activity. It integrates genetic, transcriptomic and proteomic variations with external environmental influences. However, the inherent complexity of metabolomic datasets, characterised by high-dimensional chemical diversity, wide dynamic ranges, and significant matrix effects, necessitates a rigorous suite of chemometric and bioinformatic workflows. For researchers uninitiated in computational biology, the multi-stage requirement for raw data pre-processing, signal deconvolution, and multivariate statistical modelling (such as PCA or PLS-DA) presents a substantial barrier to entry. Navigating these convoluted data architectures remains a primary challenge in deriving biological meaning from the global metabolic profile. Here, we present a workflow to use Python Dash Apps to create a user-friendly interface for simplifying data processing and statistical calculations. Users can select their desired samples to initiate calculations for various statistical tests, generating interactive and publication-quality figures to explore their results. These apps were deployed on an Apache server via cPanel, allowing individuals to share their findings with collaborators and for research facilities to share metabolomics results with their users.

Computer-aided drug design for conditional biomolecular interfaces requires evaluation across more than one receptor structure, docking pose, or scalar score. LINE-1 ORF1p is treated here as a state-family interface target whose relevant behavior is distributed across receptor microstates, assembly-compatible contact neighborhoods, ligand conformers, and perturbation snapshots. This article presents Linobectide as a mathematical-chemistry CADD workflow centered on a modified black-hole algorithm (MBHA) for persistence-weighted prioritization of putative ORF1p inhibitor candidates. Each molecule is represented as a dossier containing standardized descriptors, docking annotations, interaction-class persistence vectors, finite-action stability traces, graph-localization summaries, SPECTRAL-SAR applicability-domain records, and rank-shift diagnostics. The revised analysis emphasizes numerical reporting endpoints: fixed run parameters, baseline comparators, ablation metrics, rank stability, regeneration fractions, protected-elite fractions, and reproducibility indices. Docking is used as an annotation layer rather than as a stand-alone proof of inhibition. The framework is therefore reported as a transparent computational prioritization protocol that generates testable hypotheses for future biochemical and cellular validation, not as experimental proof of ORF1p inhibition or therapeutic activity. Author summaryDrug-design workflows can become over-dependent on the best docking pose even when an interface target remains functional through alternative contact corridors. Linobectide addresses this issue by ranking candidates only after docking annotations are aggregated across receptor-state and perturbation conditions. The MBHA search promotes a candidate when interaction persistence, finite-action stability, graph localization, SPECTRAL-SAR coherence, applicability-domain support, and reproducibility checks are concordant. The revision removes unsupported claims of performance advantage and replaces them with benchmarkable endpoints that can be compared with docking-only, consensus-docking, and ablated MBHA baselines. The SI Appendix is retained as a figure atlas for state-family construction, graph-localization diagnostics, docking provenance, consensus geometry, and comparative triage.

The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.